POLYHEDRAL APPROXIMATIONS OF STRICTLY CONVEX COMPACTA 
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Abstract. We consider polyhedral approximations of strictly convex compacta in finite dimensional Eu- 
clidean spaces (such compacta are also uniformly convex). We obtain the best possible estimates for errors 
of considered approximations in the Hausdorff metric. We also obtain new estimates of an approximate 
algorithm for finding the convex hulls. 



1. Introduction 

We begin by some definitions for a finite dimensional Euclidean space (K", || • ||) over K with an inner 
product (•, •). Let B r (a) — {x <G M" | \\x — a\\ < r}. Let cl A denote the closure and mtA the interior of 
the subset A C W 1 . The diameter of the subset A C K™ is defined as diam A = sup \\x — y\\. The distance 

x,y£A 

from the point x € K™ to the set A C R™ is given by the formula g(x, A) — inf \\x — a\\. We shall denote 

a£LA 

the convex hull of a set A C W 1 by co A, the convex hull of a function / : W 1 -> K by co / (cf. [Tl 151 [TB"] ) . 
The Hausdorff distance between two subsets A,Bc M. n is defined as follows h(A, B) = 

= max | sup inf||a-6||, sup inf lla - 611 1 = inf{r > I A C B + B r (0), B c A + B r (0)}. 

[aeA bEB beB aeA J 

The supporting function of the subset A C 1™ is defined as follows 

s(p,A) = sup(p,x), Vpel" (l.l) 
xeA 

The supporting function of any set A is always lower semicontinuous, positively uniform and convex. If the 
set A is bounded then the supporting function is Lipschitz continuous [TJ [5] . 

It follows from the separation theorem that for any convex compacta A, B in M. n (cf. [HI Lemma 1.11.4]) 

h(A,B) = sup \s(p,A)-s(p,B)\. (1.2) 

IIpII=i 

A convex compactum in R™ is called strictly convex if its boundary contains no nontrivial line segments. 

Definition 1.1. (Polyak [H]). Let E be a Banach space and let a subset A C E be convex and closed. 
The modulus of convexity 8a ■ [0, diam A) — > [0, +oo) is the function defined by 



Bs ( Xl ~l X2 ) c A i Vxi,3c 2 e A : ||a5x - ^2 1| = £ 



<5 A (e) = supjj > 

Definition 1.2. (Polyak [TT]). Let E be a Banach space and let a subset A C E be convex and closed. 
If the modulus of convexity <5^(e) is strictly positive for all e € (0,diam^4), then we call the set A uniformly 
convex (with modulus 5a{'))- 

We proved in [5] that every uniformly convex set is bounded and if the Banach space E contains a 
nonsingleton uniformly convex set then it admits a uniformly convex equivalent norm. We also proved that 
the function e — ¥ Sa(s)/s is increasing (see also [5] Lemma l.e.8]), and for any uniformly convex set A there 
exists a constant C > such that Sa(s) < Ce 2 . 

The class of strictly convex compacta coincides with the class of uniformly convex compacta with moduli 
of convexity 5a (e) > for all permissible e > in the finite dimensional case (cf. [3])- 

Definition 1.3. ([HI IS])- A grid G with step A 6 (0, |) is a finite collection of unit vectors {pi} C M. n , 
i e 1,7 = {1,...,/}, such that for any vector p ^ 0, p € K", with A ^ G there exist a set of indexes 
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I p C 1,1 and numbers at > 0, i £ I p , with the property 

p=^a;Pi, Pi G G, (1.3) 

IIPi-Pj-H < A, Vijelp. (1.4) 

It is well known [TJ El [13] that for any convex closed subset A C R n we have 

A = {xeW l | (p,z) <sb,A), VpeB x (0)}. 

We shall consider external polyhedral approximation of the compact A C R" on the grid G from Definition 
11 

i = {i£l" | (p, z) < s(p, A), Vp € G}. 

From the inclusion G C Bj(0) we easily see that A C A. For an arbitrary convex compact set A C R™ we 
have ft(A,i) < 2h({0},A)A (cf. 00). If A= f| B fi (x) ^ 0, then /i(A,i) < 2i?A 2 (cf. 0(3]). Further 

we shall consider the approximation of an arbitrary strictly ^uniformly convex compact set A C M n . Our 
further goals are 

1) Estimate the error h(A, A) via the geometric properties of the set A; 

2) Suppose that we know a presupporting function f(p) of the convex compact A, i.e. the function / is 
positively uniform, continuous and co/(p) = s(p,A). Let 

A = {x G R n | (p, x) < f( P ), Vp G G}. 

Estimate the error h(A, A) via the properties of the function / and geometric properties of the set A. In 
this case we do not know the supporting function of the set A, but we can find information about some 
properties of the set A: diameter, modulus of convexity, etc. 

3) We shall consider an algorithm for calculating the convex hull of a positively uniform function defined on 
the grid and discuss estimates for the errors of such algorithms. 

2. Approximation by supporting functions 

Lemma 2.1. For a given grid G (Definition with step A G (0, ^) in the representation of any vector 
p 7^ 0, jj^jj- ^ G by formulae the following estimates hold 

||p-^||<A, VjG/p, 1 > HIH > 1 - ^A 2 , 

where 

<p= — = S^aiPi, a=y^c<i, &i = —. (2.5) 



Proof. From the definitions of p and on we obtain that 

p = o-iPi, &i = i. 



Hence 



^2 aii(p-Pi) = o. 



By the triangle inequality we get 

IIP -Pill < X! &l W Pl _ PiH < A ' V -? e J p> 



HpII < E ^INI = L ( 2 - 6 ) 



i£l„ 
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The condition \\pi — pj\\ < A is equivalent to the condition (pi,Pj) — ^(|b;|| 2 + |bj|| 2 — \\Pi ~ Pj\\ 2 ) > 
l-A 2 /2. Thus 1 > II^H > ||p|| 2 = j:,^. ^^(p^p^) > (1 - |A 2 ) cti^- = 1 - |A 2 . " ' " □ 

Definition 2.1. Let / : R" — > M be a positively uniform function. Let G be a grid with step As (0, |). 
Define the grid operators 

( f(p), ieG, f /(?)> m e G ' 

c/(f)= W/(p) = 

i +oo, A £ G, E nfir ^ 

v iipii V ieip 

Indices I p and numbers on are from Definition 1 (3j 

Lemma 2.2. (^[8j Lemma 5], [9l Lemma 2.6.2]/ Let f : M. n -}Rka positively uniform function. 
(1) If the function f is convex then Cf > /, coCf{p) = f{p) \/p G G. 
(S) If the function f is convex then f < colAf. 

(3) coCf = coUf, V/. 

(4) co/<coC/ ; V/. 

The next lemma is a modification of Lemma 2.2 from [3]. 

Lemma 2.3. Let A C M n 6e compact and uniformly convex set with the modulus of convexity 8. Let 
e G (0,diamA), A G (0, |). Let p Xl p 2 G R", ||pi|| = 1, 1 - ±A 2 < ||p 2 || < 1. Let x 4 = argmax(pj, x), 



i=l,2. //|bi-p 2 || <(4-A 2 )%l, tften||a;i-a:2|| < e. 



Proof. Suppose that ||a;i — x%\\ > e. Let t = S(\\xi — x 2 \\) > S(e). By the condition 

c A 

we have that (pi,x±) > (p±, Xl + X2 } + t, 

(pi,x 1 -x 2 )>2t, (2.7) 
Xl+X 2 \ ^ f x% + x 2 \ , ,(, 1 a2 



(P2,x 2 ) > (p 2 , ^f^j + t\\ P2 \\ > (p 2 , ^f^j + t { 1 -2 A 

(P2,x 2 - Xl ) > (2- A 2 )t. (2.8) 

By formulae (2J7}, (2|5j! we obtain that 

(pi -i>2,aci -x 2 ) > (4- A 2 )t, 

and 

llPi-pall >(4-A 2 ): 



2 c5(||a;i — ^2 1|) 



xi - x 2 \ 

By Lemma 2.1 of [3] we have the inequality ^J^^j^ > ^— and 



|b 1 -p 2 ||>(4-A 2 )^. 

Contradiction. □ 
Corollary 2.1. Let under the conditions of Lemma 2\^e(A) be a solution of the equation — 4 _^ A -i ■ 

If \\pi - P2W < A then \\xi - x 2 \\ < e(A). 

Proof. The proof follows from Lemma 2 13] and strict monotonicity of the function [H Lemma 1.2]. □ 

Theorem 2.1. Let Pl ,p 2 G R n , \\pi\\ = 1, 1 - ±A 2 < ||p 2 || < L Ibi - Pz\\ < A. Let A C R" be 
compact and uniformly convex set with modulus of convexity 5 and A G (0, |), <5(diam A) /diam A > . 
Let Xi = argmax(pi, x), i — 1,2. Then 

s(pi,A) - s(p 2 ,A) = (x 2) pi -p 2 ) +£i(|bi -palDIbi -Pall, (2-9) 
s(p 2 ,A) - s(pi,A) = {xi,p 2 -pi) +e 2 (|bi -PalDIbi -Pali, (2-10) 

and 

max{| £l (|bi -p 2 ||)||bi -P2II, kidbi -P2||)||bi -P2II} < e(A)A, 
where e(A) is a solution of the equation ^2 = 4 _^^ . 
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Proof. Equations (2^ and f2 fTU]) are equivalent to the condition of continuous gradient for (convex) 
supporting function at the points p 2 and p\ , respectively. It is a well known fact that the supporting function 
of the strictly convex compact is continuously differentiable ( [HI H3] ) - 

By Corollary 2Q]we have the estimate 

|(pi ~P2,xi - x 2 )\ < \\pi -pallllari - x 2 \\ < Ae(A). 
From the equalities (j>i, Xj) = s(pi, A), i = 1,2, we conclude that 

EiGbi -P2IDHP1 -P2II = (Pi,xi - x 2 ), eadbi -PalDIbi -Pa || = (.P2,x 2 - x x ). 
By the formulae 

(pi,xi - x 2 ) + (p 2 ,x 2 - xi) = (pi -p 2 ,xi - x 2 ), (2-11) 
and (pi, ii — X2) > 0, (p 2 ,x 2 — x\) > we have 

max{|(pi,xi - x 2 % \{p2,x 2 - x\)\} < (pi -p 2l x 1 - x 2 ) < e(A)A. 

□ 

It is well known that the external polyhedral approximation A of the convex compact set A C M. n with 
supporting function s(p, A) on the grid G satisfies the formula s(p, A) = coCs(p, A) ([8j [9]). 

Theorem 2.2. Let A C M™ be a convex compact set with the modulus of convexity 5(e), e € [0,diamA|. 
Let & be a grid with the step A £ (0, <5(diam A)/di&mA > ■ Then 

h(A, A) < ^£(A)A, 

where e(A) is a solution of the equation = 4^37 ■ 

Proof. From the inclusion A C A, formula s(p, A) = coCs(p, A) = coUs(p, A) (see Lemma 2J5J) and from 
Definition 2 JT| it follows that (in terms of Definition 1 J3j) 

< s(p,A) - s{p, A) <Us(p,A) - sip, A) = aJ2 a t (s( Pl , A) - sip, A)), Vp e M™. (2.12) 

iei P 

Let i = argmaxfp, x). 

x<EA 

sipi,A) - sip, A) = (x,pi-p) +£i(|bi-p||)|bi-p||, 
and by properties of vector p (Lemma 2fT]) and Theorem 2Q] we conclude that 

|£i(lbi-p||)lbi-p||| <e(A)A, Vie I p . 

Finally, using the equality ^iPi — P> we obtain that 

ieip 

&i (sipi,A) - sip, A)) = Y & i i(x,Pi-p) +£i(lb» -PlDIbi -Pll) < 
< max|ei(|bt -p||)||Pi -p||| < e(A)A. 

By formula (2ll2"]) it follows 

< sip, A) - sip, A) < ae(A)A = |p|e(A)A < -^^||p|| < |e(A)A||p||. 

Ibll 1 - 5 A 7 

By formula (ill]) we get 

h(A,A)<*e(A)A. 

□ 

Corollary 2.2. For any convex compact set A the modulus of convexity 5 satisfies the estimate <5(e) < 
Ce 2 . Thus the typical value of e(A) is e(A) x A s , A — » +0, where s e (0, 1]. 
Corollary 2.3. The estimate of Theorem 2\^is exact. 

Consider an example. Let M 2 be the Euclidean plane with the standard basis Ox\x 2 . Let A — {x 2 > 
\x 1 \'}nB 1 (0),s>2. 
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The modulus of convexity for the set A equals 8(e) — e s /2 s for small e > and it is realized on the 
segment [(— §, §?), (§, §?)] ■ Let a and b be two points from A: 

2' 2?) ' ~ V2'2^ 

Let p a and pb be unit normals to the set A at the points a and b, respectively. It is easy to calculate that 
(s(s/2y-\-l) _ (-s(e/2y-\-l) 2 S ( £ /2)^ 1 



Pb = r — - == =; =, Pa = r — =; =, \\Pb-Pa\\ = 



v / s 2 (e/2) 2 (^ 1 ) + 1 V^PF^+l A /s 2 ( £ /2) 2 ( s - 1 ) + 1 

Let A = ||p6 — poll an d Pa, Pb be adjacent vectors of some grid G with the step A. Suppose that the grid G 
has symmetry with respect to the line Ox2- 

Then e x A 7 ^ (for small A). The tangent line to the graph X2 = \%i\ s at the point b is 

Ifcmfo) = »(0' (a*" |) + (|)'- 

We have yt O n(0) = — (s — 1) (f) S - The point c = (0,ytan(0)) belongs to the set A (because approximation A 
has symmetry with respect to the line Ox2)- Hence 

h(A,A) > g(c,A) = \y tan (0)\ = (s - 1) (|)\ 
So we have h(A,A) > C ■ A^ 1 . The same order is given by Theorem 2j2] 

3. Approximation by presupporting functions 

Suppose that we know a presupporting function / of a convex set A c K". We want to estimate the 
difference coCf(p) — cof(p) for all p e M™. The last question is equivalent to the question of evaluation of 
the value h(A,A) where A = {x | (p,i) < co/(p), Vp G M"}, I = {a; | < coC/(p), Vp G M"}. 

The geometric difference of sets B,ic M ra is the set 

B-A={i|s + AcB} = P| (.B — a). 

We shall obtain the solution for two particular cases of presupporting function. The first case is f(p) — 
s(p, B) — s(p, A), where A and B are convex compacta. In the case B — A^% the convex hull co f(p) equals 
the supporting function s(p, B — A) of the geometric difference B — A [9l Formula (1.11.18)]. 

The second case is f(p) = mm{s(p, A), s(p, B)}, where A and B are convex compacta. In the case 
A n B 7^ the convex hull co f(p) equals the supporting function s(p,AD B) of the intersection A n B [9j 
Formula (1.11.17)]. 

The considered cases have important role for computational geometry [5] and for linear differential games 

[una. 

Theorem 3.1. Let A, B C K™ be convex compacta and suppose that B is uniformly convex with modulus 
8 B - Let f(p) = s(p,B) - s(p,A). Let B ro (a) C B — A C B d (a). Let G be a grid with step A G (0, |); 
<5 B (diamB)/diami? > 4 ^ 2 . Then 

OJ 

cof(p)<coCf(p)<cof(p) + —e B (A)A\\p\\, Vp e M", (3.13) 

7r 

where e — Sb(A) is a solution of = 4 A A i ■ 
Proof. The left inequality in ^3 JT5|) is obvious. By Formula (1.11.18)] we have 

co f(p) = co (s(p, B) - s(p, A)) = s(p, B — A). 

Let C = B — A. By the formula 

Cf(p)+Cs(p,A)=C S (p,B) 

we obtain that 

co (Cf(p) +Cs(p, A)) = coCs(p,B) = s(p,B), 

where B — {x \ (p, x) < Cs(p, B), \/p G M™}. Using inequality co (/ + g) > co / + co g (which is true for any 
functions /, g) we have 

coCf( P )+coCs(p,A) < s(p,B), 
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or 

coCf(p) < s(p, B) - coCs( P , A) = s(p, B) - s(p, A) < s(p, B) - s(p, A). 
By the last inequality 

coCf( P ) < co (s(p, B) - s(p, A)) = s(p, B - A). 
Let \\p\\ = 1. Using (1 |2j we have 

coC/O)- co f(p) < s(p,B — A)-s{p,B — A) < h(B — A,B — A). 

Let h = h(B, B). Using conditions of the theorem we conclude that 

B^AdB^Ad{B + B h (0)) ^Ac(B+ A(C - a)) — A = 
' B + ±(B^A)) i A - ±a c (b + ±B -i — A] ^A-±a = 



= (B^)+i((fli4)-a)c(Bi4) + AB# 

Hence 

h — AS — A) < — B). 

V / r 

Applying Theorem 2j5]we finish the proof. □ 

Theorem3.2. Let A,B C R™ be uniformly convex compacta with moduli 5 a, 5b ■ Let f{p) = min{s(j>, A), s(j>, B)}. 
Let B ro (a) c AC\B, maxjdiam A, diamB} < d. Let G be a grid with step A £ (0, \); 5a (diam A)/diam A > 
4^r, 5 B (diamB)/diamB > j^i- Thencof(p)< 

< coCf( P ) < co f(p) + j Lax{ £A (A), £s (A)} + ^ (e A (A) + £ B (A))J A||p||, Vp e R", (3.14) 

w/iere e = £jc(A) is a solution of Sx ^ = 4 A A2 , X = A or X = B. 

Proof. Let C = Ap\ B and C be external polyhedral approximation of the set C on the grid G, By the 
inclusions C C A, C £ B wc have C C A, C C B and thus CciflB. 

C/(p) = mm{Cs(p,A),Cs(p,B)}, coC/(p) = co min{Cs(p, A),Cs(p,S)}. 
By the inequality co min{/, g} < minjco /, co g} (which is valid for any functions /, g) we have 

coCf(p) < min{coCs(p, A), coCs(p, B)} = min{s(p, A), s(p, B)}. 

Hence 

coC/(p) < co (min{s(p, J[),s(p,S)}^ = s|p,i(lBj. 
Let ||p|| = 1. Using (lj2j we have 

co C/(p) - co f(p) < s (p,Ar\Bj - s(p,Ar\B) < h(AtlB,AnB). 
Applying Theorem 3.1 [2] we obtain that 

h(A n B, A n B) < max{h(A, A), h(B, B)} + (h(A, A) + h(B, . 
Using of Theorem 2j5]ends the proof. □ 

4. On finding the convex hulls 

Theorem 4.1. Let A C R n be a uniformly convex compact set with modulus of convexity 5. Let 

r = sup{r > | 3 a £ R n : B r (a) C A}. 
Let a point a £ R" be such that B ro (a) C A and d = sup \\x — a\\. Then 



d<max{2r ,r + ( r 1 ) } 



where S^ 1 is the inverse function for modulus 5. 

Proof. Suppose that d > 2r . We shall prove that d < r + 5~ x (^-)- 

Let b £ A and \\a — b\\ = d. Let L be any 2-dimensional affine plane which contains points a, b. Let 
L = (a — b) , dim L = n — 1. 
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Our further consideration will take place on the plane L. Let the line I be orthogonal to the line aff {a, b} 
and a € I. Let {x,y} — I l~l B ro (a). From the triangle xab we have ||x — b\\ > d — ro, from the triangle yab 
we have \\y — b\\ > d — ro. 

Let z = and let the line l\ be parallel to the line / and z G h. Let x\ = ^naff {x, b}, y\ = /iflaff {y, b}. 
By the uniform convexity of the set A we obtain that 



z, z + 



Xi 



\X\ - z\ 



ro 
2 



hence 



z,z + 



yi - z 

- z \\ 



ro 
2 



<5(||x-6|| 



S(d-r ) 



U 



U 



z, z + 



z,z + 



yi 



yi - z \ 



X\ 



\x\ - z\\ 



2 

111 
2 



5{\\v-b\\) 



S(d-r Q ) 



c A 



C A. 



If R = 7 -j + 5(d — ro) > ro then (due to the previous inclusion being valid for any 2-dimensional plane L 
with {a, b} C L) we have 

B R (z) n(h + z) cA 

and thus 

co (B ro (a) U (B R (z) n (L + z))) c A. 

By the last inclusion and by the inequality \\a — z\\ > ro we obtain that a shift of the ball B ro (a) on a small 
distance in the direction b — a occurs in the interior of the set A. Hence ro is not the maximal radius of balls 
from A. This contradiction shows that ^ + S(d — r ) < r . □ 

Now we describe an algorithm for finding the convex hull of a positively uniform function [9l [10] . 

Suppose that G is a grid with step A S (0, f(p) is a positively uniform continuous function and 
A = {x | (p, x) < Cf(p), Vp}. We wish to calculate coCf(p) for all peG. In other words, we wish to find 
CcoCs(p 1 A). The problem can be solved as a collection of problems of linear programming: for all q £ G to 
find 

{q,x) —} max (p, x) < Cf(p), \fp € G. 

We shall describe an approximate algorithm from [7], [5J Theorem 2.6.3], |10] and discuss its error for the 
case of uniformly convex set A = {x \ (p,x) < cof(p), Vp £ 1"} with modulus of convexity 8. 
Suppose that B ro (a) C A is the ball of maximum radius in the set A and d — sup \\x — a\\. 

We often do not know the precise values of a, ro, but we can easily calculate the ball of maximum radius 
BrQ)), ro < R, from A: it suffices to solve the following problem of linear programming 

R ->■ max (p, b) + R< f{p), Mp e G. 
The solution (b, R) € M" x M gives the center of the ball and its radius. In this case Bn(b) C A and 



Acb- 



A) Bi(0), see [9j Corollary 2.6.2]. Thus 



(p,b) + R\\p\\<f(p)<( P ,b)+[d + 



4d 2 



-A 



(4.15) 



We shall further assume that 6 = 0. 

The first step of the approximate algorithm is to calculate for all q £ G the values 

s (q) =max|-^-y | p£ 

The second step is to define z(q) = q/s°(q), \/q G G. Then the polyhedron 

A\ = co (J z(q) 



is 



is an approximation of A. The approximate value for coCf(p), p G 

ma,x{(p,z(q)) | q G G}. 
By [SJ Theorem 2.6.3] we have under the assumption (4 lT5t that Ai C A, 



h(A,At) < 



(d+^A)' 



R 



A. 



(4.16) 
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and 



< coCf{p) - max{(p, z{q)) \ q G G} < 



d + —A 

r 



R 



Vp G G. 



So the error of the algorithm is proportional to the step A and to the value — in the general case. 
Consider the case when the set A has modulus of convexity of the second order: 5(e) = Ce 2 + o(e 2 ), 
-> +0. Then under assumption A G (0, y^ro) we obtain by Theorem 4Q]that 



d+^A 

ra 



< 



d+^A 

ra 



< Const. 



R r 

and the error of the algorithm does not depend on the radius R of an interior ball. 

5. Epilogue 

1. By Theorem 4[T]we can estimate the value in the theorems from Section 3. For example, in Theorem 
3Q]in the case when B To (a) is the ball of maximum radius from B — A and d — sup ||x — a\\, we have 



(for small ro > 0) 



d_ ro_ 



'ra ' 
v 2 , 



If the modulus 5g has the second order at zero then — 



\/ra 



ro 



+0. 



2. We want to point out that if the sets A, B are uniformly convex with moduli 5a, 5b respectively and 
A, B are polyhedral approximations of A and Bona grid G with the step A, then 

h(A + B,A + B) < 5 £A+S (A)A 

and in general in spaces of 3 or more dimensions A + B C A + B, but A + B ^= A + B. So the sum of 
approximations does not equal the approximation of sum. 

3. The results can easily be reformulated in any finite dimensional Banach space. The only obstacle for 
the proofs is in Lemma 2Q]when we estimate ||p|| > 1 — \A 2 . One must demand from the space and the 
grid that 



C= inf 

IHI=i 



G(0,1). 



^2 &ipi 

ieip 

Then one must replace denominator 4 — A 2 by the new 2 + 2C and coefficient | by the new ^ in all theorems. 
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